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Discrete Klein-Gordon models with static kinks free of the Peierls-Nabarro potential 
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For the nonlinear Klein-Gordon type models, we describe a general method of discretization in 
which the static kink can be placed anywhere with respect to the lattice. These discrete models are 
therefore free of the static Peierls-Nabarro potential. Previously reported models of this type are 
shown to belong to a wider class of models derived by means of the proposed method. A relevant 
physical consequence of our findings is the existence of a wide class of discrete Klein-Gordon models 
where slow kinks practically do not experience the action of the Peierls-Nabarro potential. Such 
kinks are not trapped by the lattice and they can be accelerated by even weak external fields. 

PACS numbers: 05.45.-a, 05.45.Yv, 63.20.-e 



I. INTRODUCTION 

Discrete solitons and more specifically kink-like topo- 
logical excitations are ubiquitous structures that arise in 
numerous physical applications ranging from dislocations 
or ferroelectric domain walls in solids, to bubbles in DNA, 
or magnetic chains and Josephson junctions, among oth- 
ers (see, e.g., pj for a recent exposition of relevant appli- 
cations) . The mobility of such lattice kinks is one of the 
key issues in many of these applications, especially since 
the pioneering works of [3, Q which illustrated that the 
kinematics on the lattice is dramatically different from 
the continuum analog of such equations where constant 
speed propagation is typical. Instead, on the discrete 
substrate, kinks need to overcome the, so-called, Peierls- 
Nabarro potential (PNp), constantly radiating their en- 
ergy and being eventually trapped by the lattice. The 
static PNp refers to the energy difference between a sta- 
ble inter-site centered discrete kink and an unstable, on- 
site centered discrete kink. Clearly, as a kink is travelling 
from one site to the next, it "wobbles" over this potential 
energy landscape Q. However, even though clearly trav- 
elling is intimately connected with overcoming the static 
PNp without "radiating" energy Q, this connection is 
relatively subtle and the inter-dependence of these two 
features (static PNp and travelling) still remains elusive 
0. Typically, discrete kinks traveling with finite veloc- 
ity have only been obtained for a discrete set of veloc- 
ities which makes the motion unstable with respect 
to perturbations. There exists a class of more exotic ex- 
act solutions (the so-called "nanoptera" ) where the kink 
propagates together with a plane wave having the same 
velocity @- 

While the travelling problem is extremely interesting 
in its own right, in the present work, we will start by 
examining the construction of discrete models with PNp- 
free kinks, using a simplified (quasi)static approach. Two 
classes of discrete models where static kink can be placed 
anywhere with respect to the lattice have been previously 
derived: one conserving energy |8( and another one con- 
serving momentum 0. In both cases the static kink so- 
lution can be obtained from a two-body nonlinear map. 



In the present paper we demonstrate that, in general, a 
discrete version of the first integral of the static contin- 
uum Klein-Gordon field plays the role of this nonlinear 
map. Thus we derive a wide class of such models in- 
cluding the two above-mentioned classes as special cases. 
The advantage of this approach is that the kinks are no 
longer (typically) trapped by the lattice. Instead they 
can be accelerated by even weak external fields. However, 
a note of caution should be added here. While one might 
naively expect that such solutions would be intimately 
connected with slow travelling, it has been demonstrated 
numerically that travelling solutions (when they can be 
found as e.g. in 0,03 for Klein-Gordon lattices, us- 
ing the methods of |ll|) have a sharp lower bound in 
their wave speed jl^. The existence of such a threshold 
illustrates the fact that one should be particularly care- 
ful in trying to infer features of the travelling problem 
from such "static" considerations. On the other hand, as 
the recent work of Barashenkov, Oxtoby and Pelinovsky 
demonstrates , discretizations without PNp are much 
more natural candidates for possessing travelling solu- 
tions for a isolated wave speeds (not close to zero). 

The presentation of our results will be structured as fol- 
lows. Section II will contain the setup and notations used 
for the Klein-Gordon models. Section III will present 
the general methodology for obtaining static PNp-free 
discretizations. Section IV will illustrate the connection 
to previously reported models. Section V will focus on 
the special case example of the </> 4 model, for which our 
numerical observations will be presented in section VI. 
Finally in section VII, we will summarize our findings 
and present our conclusions. 



II. SETUP 



We consider the Lagrangian of the Klein-Gordon field, 



L = 



dx . 



(1) 
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and the corresponding equation of motion, 

<kt = <f> xx ~ V'{4>) = D(x) . (2) 

Topological solitons (kinks) are possible only if V((j>) 
has at least two minima 0qi an d 002, where V'(^>oi) = 
and V"((poi) > 0. Obviously, 4> — 4>oi an d 4> = 4>02 
are the stationary solutions to Eq. J2J. We will study 
the properties of kinks that interpolate between these 
two stationary solutions. Our considerations allow one 
to treat the cases when other minima appear in between 
the two minima, 4>oi and (f>02, connected by the kink. 

Equation J2J will be discretized on the lattice x — nh, 
where n = 0, ±1, ±2..., and h is the lattice spacing. 

For brevity, when possible, we will use the notations 

<j>n-i =1, <pn = m, <j) n +i = r . (3) 

We would like to construct a nearest-neighbor discrete 
analog to Eq. J5J of the form 

m = D(C,l,m,r), (4) 

where C > is a parameter related to the lattice spacing 
h as C = l//i 2 , such that in the continuum limit (C — ► 
oo), D(C, I, to, r) D{x) = (j) xx - V'{<t>). 

Note that in this context, the "standard" discretization 
emerges in the form: D(C,l,m,r) = C(l — 2m + r) — 
V'{m). Generalizations of this model will be discussed 
in the form 

fh = C(l -2m + r) - B(l,m,r), (5) 

where B(l, m, r) has V'((j)) as the continuum limit. 

We will characterize a model as PNp-free if a static 
kink can be placed anywhere with respect to the lattice 
(continuum, rather than discrete, set of equilibrium so- 
lutions). This is equivalent to demanding that the kink 
must have an neutral direction, or (from Noether's the- 
ory 0]) a Goldstone translational mode. It is natural to 
categorize this definition of PNp-free model as "static" 
or "quasi-static" , in the sense that it does not involve the 
kinematic or dynamical properties of the model. 

On the other hand, one can demand the absence of 
PNp at finite kink velocities. This can be transformed to 
the demand that the discrete model supports the exact 
traveling wave solutions and this demand can be called 
"dynamic" definition; see e.g. for such travelling wave 
examples, where the "static" definition of the PNp clearly 
fails. 

In this paper we aim to construct the models PNp-free 
in the static sense as a first (yet nontrivial) step towards 
understanding the nature of the discrete travelling prob- 
lem (see also the comments above). 

We will also focus on the existence of physically 
motivated conserved quantities for the derived models. 
Hamiltonian models are energy-conserving models and 
the models with dM/dt — 0, where 

oo 

M = (<?Wl - <t>n-l) , (6) 



will be called momentum-conserving models. As was 
shown in |jg, the discrete model of Eq. Q conserves 
the momentum of Eq. ffif . if it can be presented in the 
form 

H(m,r) -H{l,m) 
m = — . (7) 

This can be verified by calculating 

n 

= ^Wn, <t>n+l) ~ 0n)] - 0, (8) 

a 

where we have used the fact that the terms 0n(0n+i 
■ cancel out due to telescopic summation. 

III. STATIC PNP-FREE DISCRETIZATION 

Our aim here will be to discretize Eq. J5J in a sym- 
metric way, so that the static kink solution can be found 
from a reduced first-order difference equation. Accord- 
ing to |8j, if we achieve that, then we are going to have 
a one-parameter family of solutions with the possibility 
to place equilibrium kinks anywhere with respect to the 
lattice (and hence, PNp-free in the static sense). 

The first integral of the steady state problem in Eq. 
©, (j> x — y/2V((j)) — (with zero integration constant), 
can be written in the form 

w{x) = g{ ( j> x )-g(yW{^))=Q, (9) 

where g is a continuous function. 
Our plan will then be the following: 

• discretize the first-order differential equation of Eq. 
iJSJ using a first order difference scheme w(l, to) = 
0. 

• Then express the right-hand side of Eq. @ as a 
sum of terms containing derivatives, e.g., dw/dx, 
dw/dcj), etc. 

• As a result, discretizations of such terms, e.g., 
dw/dx ~ \/C[w(m,r) — w(l,m)], vanish for 
w(l,m) = (or otherwise stated: the construction 
of the equilibrium solution is converted to a first 
order difference problem). Then, the static kink 
(PNp-free, by construction) solutions for the ob- 
tained discrete model can be found from this two- 
site problem. 

In the following, we will consider a particular case of 
Eq. 10 with g(£) = £ 2 , for which we introduce the nota- 
tion 

u{x) = <g-2V{^) = 0, (10) 
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and the following two-site discrete analog 

u(l,m) = C(m- I) 2 - 2V(l,m) = 0. (11) 
We will also use the shorthand notations, 

ui=u(l,m), u m — u(m,r). (12) 

We have assumed that the Klein-Gordon field supports 
kink solutions. Then, at least for the case of weak dis- 
creteness, Eq. (ITTT) also supports static kinks because 
it is nothing but a discretization of the first integral of 
static version of Eq. (see also @). 

The next step is then to find a discretization of the 
right-hand side of Eq. J5J, D(x), which vanishes when 
Eq. Q} is fulfilled. 

One simple possibility comes from the following finite 
difference 



Di(l,m,r) 



ui 



I 



1 du 
2dd> 



D(x) 



(13) 



One can also consider, more generally, continuous func- 
tions q(ui, h) such that q(0, h) = and, in the continuum 
limit, q(u, 0) = u and j^(u, 0) = 1. For example, one can 
take q = (e hu - l)/h or q = u + Y, n> i Aih n ~ 1 u n with 
constant A n , etc. Then, 



Then, by construction, any structure derived from the 
two-site problem of Eq. I|llfl is a static solution of Eq. 
(TT%1) and hence, the latter is the static PNp-free discrete 
model. 

The model of Eq. I|18[) can be generalized in a num- 
ber of ways. For example, function D3, Eq. (|17f) . can 
be modified choosing different functions V(l, to, r) to dis- 
cretize V(<f>). Then, the modified D3 can be added to the 
linear combination in the right-hand side of Eq. I|18l) . 

The model of Eq. I|18[) can be also generalized by ap- 
pending terms which disappear in the continuum limit 
and ones that vanish upon substituting ui = and 
u m = 0. For example, the derivative df(u)/d<p can be 
discrctized as 2[f(u m ) — ,f(ui)]/(r — I) or as 2f'(ui/2 + 
u m /2)(u rn — ui)/(r—l). If then we have difference of such 
terms in the equation of motion, then in the continuum 
limit they will cancel out. 

Any term in the right-hand side of Eq. I|18|l can be 
further modified by multiplying by a continuous function 
e(C, I, to, r), whose continuum limit is unity (see e.g. 
for such an example, also discussed in more detail below). 

Generally speaking, the discrete PNp-free Klein- 
Gordon models derived here do not conserve either an 
energy, or a momentum-like quantity. However, as it will 
be demonstrated below, they contain energy-conserving 
and momentum-conserving subclasses. 



1 dq ( dq 

2 d(h \du 



D(x). 



(14) 



Discretizing the left-hand side of Eq. Ijl4|l we obtain 

1 1 



lq(u mi h) - q(u u h) 
D2 =2 — 



(15) 



_q'{ui) q'{u m ) 
Inspired by [8j, we note that, in the continuum limit, 
v(m,r) v(l,m) dv <p xx 



and find 



£3: 
'y/C(r 



I 



a, 



(16) 



-I 



y/2V(l,m,r) x 



III 



TO — I 



du 

d~4> 



'2V 



dV2V 



<2V- 



D{x). (17) 



Since the expressions for Di(l,m,r) given by Eqs. 
(|13fl . (|15|) and 117|) tend to D{x) in the continuum limit, 
one can write the following discrete analog to the Klein- 
Gordon equation Eq. @ 

to = y ] bjDj(l, m, r), where Vjfej = l. (18) 



IV. CONNECTION WITH PREVIOUSLY 
REPORTED MODELS 

One energy-conserving PNp-free Klein-Gordon model 
has been derived by Speight with co-workers [8j with the 
use of the Bogomol'nyi argument [l(j. Their model, can 
be written in the form of Eq. JSJ , with the Lagrangian 



L 



\ ^ 12 



C 



2 ^ r n 



<l>n-l) 



1 1 



y-y ( G{4> n ) - G((f> n -i) 



n 



where G'(cj>) = y/V {<)>). 



(19) 



The static kink solution can then be derived from the 
lattice Bogomol'nyi equation 8], which can be taken in 
the form 

U(l, m ) = C(m-ir-2( G{m) -f il) ) 2 = 0, (20) 

\ TO — I J 

which is a particular case of Eq. 111|) . The equation of 
motion derived from Eq. (|19|) . written in terms of Eq. 
(HU, is 



U„ 



Ui 



r — to to — I 



y / 2V(m) x 



^C(r - to) 2 - U m ^C{m - l) 2 - U t 



(21) 
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The right-hand side of Eq. (|21|l is a particular case of 
D 3 (l, to, r) given by Eq. JTJJ with V(l, m, r) = V(m). 

Momentum-conserving PNp-free models were pro- 
posed in and further studied in are the 
non-Hamiltonian models of the form 



to = Di(l, to, r), 



(22) 



where D\ is given by Eq. (fT^|) . Notice that Eq. (|2~2")l can 

be mapped into the formulation of Eq. . Static kink 
solutions in this model can be found from Eq. (|1 1|) . 

If Eq. (fill) is taken in the particular form of Eq. (|20|) . 
then the momentum-conserving PNp-free model Eq. 122|) 
and the energy-conserving PNp-free model Eq. I|21|) have 
exactly the same static kink solutions. It has been proved 
that a standard nearest-neighbor discrete Klein-Gordon 
model conserving both energy and momentum does not 
exist E3. 



APPLICATION TO THE 



MODEL 



As an example, we will discretize the well-known </> 4 
field theory with the potential 



V(cj>) 



(23) 



By construction, the PNp-free models derived above 
are written in singular form. In this form the equations 
are inconvenient in practical simulations and one may 
wish to find such particular cases when singularities dis- 
appear. 

For example, for the energy-conserving PNp-free model 
expressed by Eqs. (fH3)) - lt2"T|) . singularity always disap- 
pears when G{4>) is polynomial Q. Particularly, for the 
</> 4 model with the potential Eq. (|23[l . one obtains from 
Eq. (|21|l the following_energy-conserving PNp-free dis- 
cretization derived in 



C + - \ (l + r - 2m) + to 



-— [2m 3 + (to + Z) 3 + (to + r) 3 ] , 
18 

whose static kink solution can be found from Eq. 
which, for the 4> 4 potential, obtains the form 

■I) + to 2 + lm + I 2 - 3 = 0. 

Now let us turn to the momentum-conserving model. 
Substituting Eq. JTJ} into Eq. we obtain 



(24) 

(29, 

(25) 



C(r - 2m + I) 



V(m, r) — V(l, to) 



(26) 



To remove the singularity, V(l, to) should be taken in the 
symmetric form V(l,m) = V(m, I), e.g., as 

V(l, to) = (1/4) - (a/2)(m 2 + I 2 ) + (a - l/2)mZ 

+(/3/2) (to 3 + I 3 ) - (p/2)ml (to + I) 

+(7/2) (to 4 + Z 4 ) + (S/2)ml (to 2 + I 2 ) 

-(j + 6- 1/4)to 2 Z 2 , (27) 



with free parameters a, [3, 7, and <5. In the continuum 
limit, when I — > to and r — > to, Eq. 127|) reduces to 
V(4>). Substituting Eq. J2IJ into Eq. d§ we obtain 
the following momentum-conserving PNp-free <fi 4 model 
derived in |17J 

to = (C + a) (I — 2m + r) + to 
-(3(l 2 + lr + r 2 ) + j3m(l + r + to) 
-7(Z 3 + r 3 + l 2 r + lr 2 ) - 5m(l 2 + to 2 + r 2 + Ir) 

+ (2j + 2S-l/2)m 2 (l + r). (28) 

The momentum-conserving model Eq. I|28|) with a = 
(3 = 7 = 5 = can be written in the form 



7/f 



2C 



C(l - 2to + r) + to - to 3 



(29) 



The following energy-conserving model, studied in [lE 



TO 



C(l -2m + r) + 



1 -to 2 /(2C)' 



(30) 



has the same continuum limit as model Eq. (|29|l . Fur- 
thermore, it can be derived from Eq. I|29|l by multiplica- 
tion with a factor e(C, l,m,r) = 1/(1 — m 2 / (2C)), which 
possesses a unit continuum limit. The model Eq. I|29() 
is PNp-free and thus, model Eq. (|30|) is also PNp-free 
since they have the same static solutions derivable from 
C(m— I) 2 — (1— ml) 2 /2 = 0. Thus, we have another exam- 
ple when energy-conserving and momentum-conserving 
PNp-free models have exactly the same static kink solu- 
tions. 

It is interesting to note that the energy-conserving 
model of Eq . 130fl cannot be constructed by the method 
reported in [8| where discretization of the anharmonic 
term always involves <\> n -\ and 4>n+\- More generally 
than it is done in Q , the problem of finding the energy- 
conserving PNp-free models can be formulated as follows. 
We need to discretize the potential energy of the La- 
grangian Eq. P in a way that the corresponding equa- 
tion of static equilibrium is satisfied when Eq. ljUJ) is sat- 
isfied. Both energy-conserving models discussed above 
are the solutions of this problem. 

As an example of model conserving neither energy, nor 
momentum we take Eq. (|15|l for the case of q(u, h) — 
u + Ahu 2 with constant A and obtain 



i m — ui (1 + Ahui + Ahu m ) 2 
r-l (1 + 2 Ahui) (1 + 2 Ahu m )' 



(31) 



This model can be obtained from the momentum- 
conserving model defined by Eq. (|13fl by multiplying 
by another function that reduces to unity in the contin- 
uum limit (/i — s- 0) . Obviously, the original momentum- 
conserving model and model Eq. (|31|l have the same 
static kink solutions. It can be demonstrated that these 
two models also have the same spectra of small ampli- 
tude vibrations and the same frequencies of kink internal 
modes. 
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VI. NUMERICS 



In our recent work |17| , some properties of kinks were 



compared for the "standard'' 
cretization having PNp, 



energy-conserving 



m = C(l + r — 2m) + m 



dis- 



(32) 



with the PNp-free models conserving energy Eq. 1(24(1 
and momentum Eq. (|29|l . 

It was found that the mobility of kinks in the PNp- 
free models is higher and also that in the momentum- 
conserving, PNp-free models a kink self-acceleration ef- 
fect may be observed. The origin of the effect is the non- 
conservative (non self-adjoint) nature of the model which, 
however, can be noticed only for asymmetric trajectories 
of particles when kink passes by [jj]]- If the trajecto- 
ries are symmetric, there is no energy exchange with the 
surroundings and kink dynamics is the same as in energy- 
conserving models, e.g., the kink self-acceleration effect 
disappears. Kinks in some of the momentum-conserving 
models was found to have internal modes with frequencies 
above the phonon spectrum. Such modes do not radiate 
and they can have large amplitudes storing a considerable 
amount of energy. 

Here we present/compare results for the energy- 
conserving PNp-free model Eq. 1(30(1 and the PNp-free 
model of Eq. 1)31(1 . generally speaking, conserving nei- 
ther energy nor momentum. For the latter model we 
take iii in the form of Eq. where the <fi 4 potential 

is discretized according to Eq. (|2*7|) and, for the sake of 
simplicity, we set a = /3 = j = S = 0. We obtain 



1 - — I C(l - 2m + r) + m - m 
Zl; J 

(1 + Ahui + Ahu m f 



(l + 2Ahui)(l + 2Ahu m ) 
where u t = C(m - I) 2 - (1 - ml) 2 /2 . (33) 

For A — 0, Eq. 1(33(1 coincides with the momentum- 
conserving model Eq. (21 • 

In the model Eq. 1(33(1 . the static kink solutions, 
phonon spectra, and frequencies of kink internal modes 
are A-independent. The energy-conserving model Eq. 
((30(1 has the same static kink solutions as model Eq. 
((33(1 but their spectra are different. The linear vibra- 
tion spectrum of the vacuum for Eq. 1(33(1 is lo 2 — 
2 + {AC - 2)sin 2 («/2) and that for Eq. CP is lo 2 = 
4C/(2C - 1) + 4C sin 2 (re/2), while the one for the classi- 



cal model Eq. ((52jl is 



2 + 4Csin 2 (K/2). 



The top panels of Fig. nP resen t the boundaries of the 
linear vibration spectrum of the vacuum (solid lines) and 
the kink internal modes (dots) as the functions of lattice 
spacing h for (a) the classical 4> A model of Eq. I(32|) , (b) 
the PNp-free model of Eq. ((30(1 conserving energy, and 
(c) the PNp-free model of Eq. I(33|) at A — conserving 
momentum. In PNp-free models kinks possess a zero fre- 
quency, Goldstone translational mode. Since all models 




5.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

t [X10 5 ] ( [x10 5 ] / [x10 5 ] 



FIG. 1: Upper panels: boundaries of the linear spectrum of 
the vacuum (solid lines) and kink internal mode frequencies 
(dots) as functions of the lattice spacing h — 1/VC. Lower 
panels: time evolution of kink velocity for different initial 
velocities and h = 0.7. The results are shown for (a) classical 
4> A model, Eq. J32J|, (b) PNp-free model conserving energy, 
Eq. (0, and (c) PNp-free model conserving momentum, Eq. 
at A = 0. 




FIG. 2: The kink velocity in the regime of steady motion [see 
bottom panel in Fig. (c)] for the PNp-free (f> 4 model Eq. 
I(3H(I is shown as a function of parameter A. For |j4| > 0.2, 
the kink self-acceleration effect disappears. 



presented in Fig. ^ share the same continuum 4 limit, 
their spectra are very close for small h(< 0.5). 

The bottom panels of Fig. ^ show the time evolution 
of kink velocity for the corresponding models at h = 0.7 
for kinks launched with different initial velocities. To 
boost the kink we used the semi-analytical solution for 
the normalized Goldstone mode, whose amplitude serves 
as a measure of the initial kink velocity. One can see 
that the mobility of kinks in the PNp-free models shown 
in (b) and (c) is higher than in the classical model having 
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PNp and shown in (a). In the energy-conserving models 
shown in (a) and (b), the kink velocity decreases mono- 
tonically due to the energy radiation. Non-Hamiltonian 
momentum-conserving model in (c) shows the effect of 
kink self-acceleration discussed in |l7| . 

It is interesting to study what happens when the pa- 
rameter A in Eq. (|33fl deviates from zero and the 
conservation law of the model (momentum conserva- 
tion) disappears. We found that the effect of kink self- 
acceleration, which can be seen in the bottom panel 
of Fig. [T] (c) for A = 0, remains for \A\ < 0.2 but 
the value of the kink velocity in the steady motion 
regime decreases with increase in \A\ as it is presented 
in Fig. |21 For \A\ > 0.2 kink self-acceleration effect 
disappears and kink velocity gradually decreases with 
time. From the above, we infer that properties such as 
the self-acceleration (for momentum-conserving models) 
or the Bogomol'nyi bounds (for energy-conserving dis- 
cretizations) render such models rather special within the 
broader class of PNp free models. However, the critical 
ingredient for the more general feature of (static) PN ab- 
sence exists in the form of a reduction of the second order 
problem into a first order. 

VII. CONCLUSIONS 

A general procedure for deriving discrete Klein-Gordon 
models whose static kinks can be placed anywhere with 
respect to the underlying lattice was described. Such 
models are called static PNp-free models. It was demon- 



strated that the models of this kind derived earlier 
H> 0, 0, ^3 are special cases of the wider family of mod- 
els derived here. 

Static kink solutions for the PNp-free models can be 
found from the nonlinear algebraic equation of the form 
u(l, m) = 0, which is a discrete analog of the first integral 
of the static continuum Klein-Gordon equation of mo- 
tion. This ensures the existence of static kink solutions 
at least for the regime of sufficiently weak discreteness 
and smooth background potential. The range of the dis- 
creteness parameter supporting stable static kinks varies 
according to the specific properties of the model. 

In this paper we have discussed only nearest-neighbor 
discretizations. However, one can easily write down a 
PNp-free model involving more distant neighbors by re- 
placing Eq. (|18fl with higher-order finite difference oper- 
ators approximating Eq. J3J, keeping the two-point ap- 
proximation, Eq. I|llfl . for the first integral of Eq. <|10l) . 

Discrete kinks in the static PNp-free models possess 
the zero-frequency translational Goldstonc mode and 
they can (almost) freely move with at least infinitesi- 
mally small velocity. Such kinks are not trapped by the 
lattice and they can be accelerated by even weak external 
fields. 

As a topic for future studies, it would be interesting 
to find any possible relation between models constructed 
here and models that support traveling kink solutions for 
finite kink velocity. Such connections are apparently un- 
der intense investigation |13| and should provide a frame- 
work for understanding travelling in dispersive lattice 
systems. 
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